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Abstract 

We consider a model of dark matter whose most prominent signature is a monochro- 
matic flux of TeV neutrinos from the galactic center. As an example of a general sce- 
nario, we consider a specific model where the dark matter is a fermion in the adjoint 
representation of a hidden SU{N) gauge group that confines at GeV energies. The 
absence of light fermionic states in the dark sector ensures stability of dark matter 
on cosmological time scales. Dark matter couples to the standard model via the neu- 
trino portal, that is, the singlet operator HL constructed from the Higgs and lepton 
doublets, which is the lowest dimensional fermionic singlet operator in the standard 
model. This coupling prompts dark matter decay where the dominant decay channel 
has one neutrino (and at least one dark glueball) in the final state. Other decay chan- 
nels with charged standard model particles involve more particles in the final state 
and are therefore suppressed by phase space. In consequence, the standard indirect 
detection signals like gamma-ray photons, antiprotons and positrons are suppressed 
with respect to the neutrino signal. This coupling via the neutrino portal is most 
robustly constrained by Super-Kamiokande, which restricts the dark matter lifetime 
to be larger than 10^^ seconds. In the near future, the scenario will be probed by the 
new generation of neutrino telescopes. ANTARES will be sensitive to a dark matter 
lifetime of order 10^^ seconds, while IceCube/DeepCore can probe a lifetime as large 
as 10^'' seconds. 



1 Introduction 



There is widespread hope that the next few years will see the detection of a dark matter 
particle. Signals of dark matter are searched for obliquely in colliders, directly in under- 
ground detectors, and indirectly in cosmic rays. Of these three methods, the last one is a 
priori the most challenging because backgrounds from astrophysical processes are difficult 
to estimate. Nevertheless, indirect detection is a topic of great current interest due to 
rapid instrumental progress and a constant flow of new and exciting data. During the past 
year, interest in indirect detection has been further amplifled by tantalizing signals from 
the PAMELA [1], ATIC [2] and FERMI [3j experiments, which may or may not be hints 
of dark matter. 

The main search channels for indirect dark matter detection are cosmic ray photons, 
antiprotons, positrons, and neutrinos. Certain features of neutrinos make them an espe- 
cially clean detection channel. Neutrinos, like photons, have no electric charge and traverse 
our galaxy along straight paths. Therefore the flux of neutrinos observed at the Earth is 
determined in a straightforward way by the production mechanism and is free of the astro- 
physical uncertainties which plague the propagation of charged particles. Neutrinos have 
a further advantage over photons in that astrophysical neutrino backgrounds are far less 
severe. Detection of a large diffuse flux of astrophysical neutrinos would thus constitute a 
smoking gun for dark matter. 

In spite of these attractive features, neutrinos are seldom singled out as a primary dark 
matter discovery channel. Neutrinos are weakly interacting and require large, dedicated 
experiments to detect. It is much easier to detect charged particles and photons. Typical 
models of dark matter encountered in the literature predict a cosmic neutrino flux that 
is smaller than or comparable to the fluxes in other observable channels. Most often, 
dark matter dominantly annihilates or decays into other standard model particles — W 
bosons, for instance, or tau leptons — which produce neutrinos when they subsequently 
decay. But then one would expect to flrst detect that dark matter signal by observing 
antiprotons from W decays or gamma ray photons from vr^s in the tau lepton decay chain. 
Annihilation or decay to primary, hard neutrinos typically appears as only one of several 
possible flnal states. The SU (2) gauge symmetry of the standard model generically implies 
that charged particles are produced at comparable rates, and those lead to more accessible 
signals. For these reasons measurement of the cosmic neutrino flux sometimes provides 
useful constraints on dark matter models [H |3] , or may serve to break a degeneracy among 
different models once a signal is detected in another channel, but is not considered the 
main discovery channel. 

In this paper we argue that one can construct dark matter models which predict a cosmic 
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neutrino flux that is more prominent than the antiproton, positron and gamma ray fluxes. 
If such a model is realized in nature, the indirect detection signal could first show up in 
neutrino telescopes. The neutrino flux from the galactic center can well be larger than the 
atmospheric neutrino background, leading to a sharp smoking-gun signature. The current 
generation of neutrino telescopes (IceCube and ANTARES) or the future one (KM3NET) 
may become the discovery machines. 

What kind of set-up is required for dark matter to produce a large neutrino flux? It 
is plausible that the dark matter particle is a perfect singlet under the standard model 
gauge interactions. If such is the case it can couple to the visible sector via gauge singlet 
operators which can be constructed out of standard model fields. The lowest dimensional 
singlet operators are the Higgs mass \H\^ and the hypercharge field strength -B^,/, both of 
mass dimension 2. Dark matter coupling through the Higgs portal or via kinetic mixing 
with hypercharge has been explored in numerous models. However, if the dark matter 
particle is a fermion, then the lowest dimensional standard model operator it can couple 
to is the dimension 5/2 operator 

HL = eabH^L^ = H+el - H^*ul = G+e" - -^{v + h + iG")uL (1) 

V2 

where v = 246 GeV is the standard model Higgs vev. We refer to this operator as the 
neutrino portal^ 

Now, suppose that the fermionic dark matter particle A is a part of a larger dark sector 
that is neutral under the standard model gauge group. Assume that the dark sector couples 
to the standard model via the neutrino portal as 

Cint = 0,ark{X){LH) (2) 

where OdarkW is a fermionic gauge singlet operator constructed from the dark matter 
particle A and other fields in the dark sector. Suppose in addition that the dark sector 
contains bosonic states that are lighter than A. If A is the lightest fermionic state in the 
dark sector, it is stable within the dark sector. But in the presence of the coupling Eq. ([2]) 
the dark matter particle is allowed to decay into fermionic states in the standard model. 
The leading dark matter decay channel is then 

A ^ (Dark) + u (3) 

which leads to a neutrino flux as the dominant signal from dark matter decay! Moreover, if 
there is only one dark sector particle in the final state, the neutrino flux is monochromatic. 



^The coupling through the neutrino portal has been previously utilized in models of dark matter in a 
somewhat different context [5]. 
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Of course, the coupling in Eq. ([2]) also predicts other, subdominant decay channels: A — 
(Dark) + {hh',W^e^ , Zu), which imply charged particle fluxes. But these other decays 
have one more particle in the final state, and are therefore suppressed by phase space. The 
branching ratio for these states with on-shell bosons depends on the multiplicity of dark 
particles in the final state and the mass of the dark matter particle, and can be made as 
small as 10~^. Cosmic neutrinos then become the most prominent signal of dark matter 
decay. 

We present here a simple model with a strongly coupled dark sector which explicitly 
realizes this idea. The model contains an SU{N) dark gauge group plus a massive Majorana 
fermion in the adjoint representation. With this low energy spectrum the dark gauge group 
is asymptotically free for any A^. Therefore the dark color is confined at low energies and 
the physical states are dark color singlets: glueballs and glueballinos (bound states of 
the adjoint fermion and gluons). The lightest, ground state glueballino (in the following 
simply called the glueballino) is our dark matter candidate. Note that the glueballino, 
being fermionic, cannot decay solely into glueballs, even if it is heavier. We assume that 
the dark confinement scale is in the (sub)-GeV ballpark which sets the mass scale for the 
glueballs, while the glueballino mass is set by the Majorana mass of the adjoint fermion of 
order 1-5 TeV. In the absence of any coupling between the dark sector and the standard 
model the glueballino would be perfectly stable. Decay channels open up once the adjoint 
fermion is coupled to the standard model via the neutrino portal, that is, via the dimension 
six coupling -r^Tr[G„^A](T„^(Li7). For A^v ~ 10^^ - 10^^ GeV and c ~ l/ievr^ (a one-loop 
factor) the dark matter lifetime is much longer than the age of the universe, while at the 
same time the decay rate is large enough to yield observable signals. The leading decay of 
dark matter is the two-body decay into a dark gluon and a neutrino. This results in a flux 
of monochromatic neutrinos from the galactic center which can be detected in the current 
generation of neutrino telescopes, ANTARES and IceCube. Other indirect detection signals 
are suppressed, so that neutrino telescopes are singled out as the most sensitive probes of 
our scenario. Moreover, the cross-section of our dark matter particle on nuclei is small, 
giving null predictions for direct detection experiments. 

The paper is organized as follows. In section 2 we introduce our example model of the 
dark sector and discuss the spectrum after the dark gauge group confines. In section 3 we 
study the thermal relic abundance, including the Sommerfeld effect and late annihilation 
of glueballinos. Depending on assumptions about modeling nonperturbative contributions 
to the annihilation cross-section, the correct relic abundance may be obtained for dark 
matter masses in the range 1-10 TeV provided the strong coupling scale is sufficiently 
small. Besides the sizable theoretical uncertainties due to modeling the strongly coupled 



3 



effects at energies below the confinement scale, the precise value of the mass is also sensitive 
to the number of degrees of freedom in the hidden sector - the larger the number of colors, 
the smaller the required dark matter mass. In section 4 we compute the decay width and 
the branching ratio for dark matter decays due to the coupling via the neutrino portal. 
We demonstrate that for the dark matter mass close to a TeV the neutrino decay channel 
is by far dominant, while other decay channels have branching ratios at a few percent 
level. For heavier dark matter, close to 5 TeV, the decay rate into charged standard model 
particles becomes comparable to the decay rate into neutrinos. In section 5 we make a slight 
detour to explain why the glueballs in the dark sector cannot be stable. We enumerate 
various options for their decay into the standard model final states and indicate options 
for model building. Glueball decays may do lead to additional indirect detection signals 
which however are highly model dependent. For this reason we choose to focus on the more 
robust signals from the glueballino decays via the neutrino portal. 

In the later sections of the paper we turn to discussing observational consequences of 
the general scenario where structure in the dark sector enables dark matter decay to yield a 
monochromatic neutrino fiux as its primary signal. In section 6 we determine the parameter 
space consistent with indirect detection experiments. The Super-Kamiokande bounds on 
the neutrino fiux from the galactic center set a lower bound on the dark matter lifetime, 
or equivalently on the scale Ajv which suppresses the coupling through the neutrino portal. 
That scale turns out to lie in the 10^^ GeV ballpark — the typical scale for lepton violating 
interactions. We also study the current bounds from the antiproton, positron and gamma 
ray fiuxes. We conclude that even for the maximum lifetime allowed by Super-Kamiokande 
the indirect bounds from other detection channels are less constraining in a large portion of 
parameter space. One exception are the bounds from antiprotons, which arc produced in 
subdominant dark matter decay channels with standard model gauge and Higgs boson in 
the intermediate state. These bounds are comparable to those from Super-Kamiokande for 
dark matter masses around 1 TeV, and when the dark matter mass is close to 5 TeV (and 
thus the branching ratio for decays into charged particles becomes large) the antiproton 
flux measured by the PAMELA experiment implies a more severe constraint on the dark 
matter decay rate than the Super-Kamiokande bounds. However, one should note that 
uncertainties involved in the propagation of antiprotons in our galaxy are very large, and 
by changing the parameters of the propagation model one can reduce the antiproton fiux by 
a factor of five, thus obliterating the antiproton bounds. We therefore continue to refer to 
the Super-Kamiokande bounds, which are comparatively free of astrophysical uncertainties. 
In section 7 we study the discovery prospects in ANTARES and IceCube. We compute 
the rate of signal events in both detectors induced by the predicted neutrino fiux from the 
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galactic center, and compare it to the background rate from atmospheric neutrinos. We 
conclude that neutrino telescopes have a great potential to probe this scenario, and their 
sensitivity will exceed that of Super-Kamiokande by 1-2 orders of magnitude. 

2 Model of Dark Sector 

In this section we construct an example of a dark matter model which leads to an enhanced 
neutrino signal. The model falls into the hidden valley class [7], where the hidden sector 
contains fairly light states (here GeV scale) interacting very weakly with the standard 
model through higher dimensional operators. The hidden sector here consists of a SU (N) 
gauge group together with a fermion in the adjoint representation. A". Incidentally, this 
is also the field content of pure Af =1 super- Yang-Mills theory, but we do not consider 
supersymmetry here: the fermion has a TeV scale mass mx and the interactions of the 
hidden sector with the standard model are not supersymmetric. 

The dark gauge group is asymptotically free for any so that the theory confines at 
low energies. We assume that the confinement scale is in the GeV ballpark, well below the 
adjoint fermion mass. The low-energy theory below a TeV contains only the dark gluons, 
and after confinement the low-energy degrees of freedom are stable glueballs. The spectrum 
of a pure Yang-Mills theory has been studied on the lattice for the case = 3 [8j. These 
studies indicate there are 12 stable glueballs in a pure SU{3) theory, labeled by their J^"-^ 
quantum numbers. The masses of the glueballs are set by the strong coupling scale A. The 
lightest glueball is the 0"*"^, with mo++ = 4.2r(^^ where Tq^ is the "force radius" related 
to the confinement scale A (defined by the running coupling in the MS-bar scheme) by 
A = 0.62ro ^ [9]. Thus, the hghtest glueball mass 

mo++ ^ 7 A (4) 

is appreciably larger than the confinement scale. The spectrum of glueball states for A^ = 3 
is shown in figure [H While lattice results for the glueball spectrum are only available for 
A^ = 3, the masses of glueballs do not vary appreciably with A^ in the large- A^ limit. The 
spectrum of the A^ = 3 theory is thus a useful benchmark for more general cases. 

In addition to glueballs, the spectrum of the dark sector contains color singlet bound 
states of the adjoint fermion and gluons. Borrowing from supersymmetric jargon we refer 
to the lightest of such states as the glueballino. The glueballino is stable under decays 
within the dark sector, and is our dark matter candidate. 
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Figure 1: Glueball spectrum for pure SU{3) gauge theory, after [8]. Masses are shown in 
units of the hghtest glueball mass ttiq. 

3 Relic Abundance of Dark Matter 



In the early universe at T > 1 TeV the adjoint fermions and the dark gluons were in 
thermal equilibrium. We assume that at even larger temperatures the dark sector was in 
thermal equilibrium with the standard model, but near a TeV the two sectors are already 
decoupled. Below a TeV the adjoint fermions become non-relativistic and drop out of 
thermal equilibrium, leaving a relic abundance that is in principle determined by the fermion 
mass, the dark gauge coupling and the number of dark colors. The freezeout calculation 
of A abundance is closely analogous to the well-studied problem of gluino dark matter 
|10[ [TT]. In the following we estimate the range of fermion masses that are compatible with 
the observed dark matter abundance. 

We first compute the relic abundance ignoring non-perturbative effects, in which case the 
standard textbook procedure applies. The perturbative matrix element for the annihilation 
A A — ^ gg can be found in |12]. At low energies the thermally averaged annihilation cross- 
section is given by 



3N^ -Ka 



where is the dark coupling at the scale rrix, given in terms of the confinement scale 
and the number of dark colors by ad = 67r/llA^log(mA/A). Integrating the Boltzmann 
equation relates the annihilation cross section to the relic abundance Yoo- The latter can 
be translated into the dark matter density fraction today by Qx ~ Yoo1^\Sq / Pmt where 
So = 2889.2/cm3, p^^t = 3.95 x IQ-^^GeV^; WMAP fixes Vtx to be 0.23 [13j. The final 
result strongly depends on the number of dark colors (the number of degrees of freedom in 
the hidden sector). Our results are plotted in Fig. [2l 
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Figure 2: Contours of relic density = -23 in the mx-A/mx plane, for = 3 (red), = 4 
(green), = 5 (cyan), = 6 (blue). Using the perturbative annihilation cross-section. 
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Nonperturbative contributions to the annihilation cross-section are significant and in- 
troduce large uncertainties [TOl [TT], [H]. Below we estimate the effects of Sommerfeld en- 
hancement and of enhanced annihilation after confinement. 

The exchange of dark gluons between the nonrelativistic fermions leads to an enhanced 
annihilation cross section at low velocity [15]. The Coulomb interaction of two massive 
particles takes place on a time scale set by the energy of the scattering particles, txx ~ 
{m\v'^)~^ ~ 1/T. Meanwhile, the time scale for a massive particle to interact with a gluon 
in the thermal bath, and therefore to randomize its color, is set by its mean free path in 
the thermal bath. This can be estimated by Txg = i = {ng{ag\v))^^ ~ m\/T'^. Since 
taa "C Txg, the two gluinos will remain in a state of definite total color throughout the 
interaction, and the sign and strength of their interaction will depend on that particular 
color state. We can solve the Schrodinger equation in the basis of definite total color, and 
derive the Sommerfeld enhancement in each color representation subspace i. Two particles 
of equal mass moving in a non-Abelian Coulomb potential have a scattering amplitude 
enhanced from the plane wave result by the factor 

E{Vr) = (6) 

where Vr is the relative velocity of the scattering particles, and ^ depends on the quadratic 
Casimir operators of the single-particle representations ri and r2 and the total representa- 
tion r, 

^ = -\{C,{r)-C,{r,)-C,{r,)). (7) 
The annihilation cross-section during freeze-out is then multiplied by the factor 

Eann{vr\ {r}) = , , . V dim{ri)Ei{vr), (8) 

dim(ri)dim(r2) ^ 

which enhances annihilation in attractive channels and exponentially suppresses annihila- 
tion in repulsive channels. 

As the gluinos are identical particles, we need to ensure that the total wave function 
is antisymmetric. Annihilation proceeds dominantly in the s-wave, where the wave func- 
tion is spatially symmetric, and therefore the net color and spin wave function must be 
antisymmetric. Thus for N = 2, the total enhancement factor E{v) is explicitly 

1 / ^"d TrOd "^"li \ 
E(Vr) = — + 9 - 5 ^-^^^^ 1 (9) 

Incorporating the Sommerfeld enhancement has an order 30% effect on the mass of the 
dark matter particle relative to perturbative freezeout, as demonstrated in figure [31 
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Figure 3: Relic density Q\ as a function of m\, with (dashed) and without (sohd) Sommer- 
feld enhancement, for N = 2. The sohd hne indicates Qx = .23. Rehc densities are shown 
for mx/A = 0.5 x 10^ (red), m^/A = 1 x 10^ (green), and rux/A = 2 x 10^ (blue). 



More difficult to estimate are the effects of the increased annihilation cross-sections after 
confinement, when the A" have hadronized [ini [HI [ll]. The heavy fermion A is localized 
within a region of order its de Broglie wavelength, l/lmxv), while the total radius of the 
glueballino is of order 1/A. While the total interaction cross-section between hadronized 
adjoints should be geometric, atot ~ ^r/A^, the annihilation cross-section should be signifi- 
cantly smaller. Estimating the annihilation cross-section is a subtle problem. The core of 
the question is whether or not two hadronized A's in a state with large angular momentum 
are kinematically able to radiate away that angular momentum and annihilate. Our hidden 
sector model differs from QCD in that there are no light fundamental quarks, and hence 
the mass gap in our hidden sector is larger compared to the strong coupling scale than 
the mass gap in QCD. We expect that this suppresses the annihilation cross-section in this 
model, relative to (e.g.) gluino freezeout, as the larger mass gap makes it more difficult to 
radiate away angular momentum. 

We therefore model the late-time annihilation cross-section by taking it to saturate the 
s-wave unitarity bound, 

(Tnp ~ (10) 

We expect that this estimate gives an upper bound on the cross-section, and therefore 
a lower bound on the relic abundance. This cross-section yields efficient annihilation of 
hadronized adjoints, lowering the relic density by roughly an order of magnitude relative 
to the perturbative result. It is therefore difficult to obtain the relic abundance Qx = 0.23 
without going to corners of parameter space, either large or m, or small A. 
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As we will explain in the following section, the most interesting choice of dark matter 
masses from the point of view of neutrino signals is the range 1-5 TeV. This range is 
consistent with the perturbative computation of the relic abundance. However, if the 
non-perturbative effects enhancing the annihilation cross-section are close to our estimate 
( JTOj) . this mass range requires either sub-GeV A or a very large (A^ > 20) number of dark 
colors. Building a complete model of the dark sector with A significantly less than a GeV 
requires some additional structure in the dark sector to be consistent with all cosmological 
constraints, as we discuss in section [51 



The low-energy lagrangian that governs the hidden sector has an accidental Z2 symmetry 
acting as A ^ — A which ensures stability of dark matter. However, general interaction 
terms do not have to respect this Z2 and may lead to dark matter decay. In this section 
we consider the coupling of dark matter to the standard model via the neutrino portal and 
study possible decay channels. 

The operator X^a^^C^^^ transforms as a fermion under the Lorentz group and as a 
singlet under both the dark and the standard model gauge groups. Therefore it has the 
right quantum numbers to couple to the standard model through the neutrino portal. The 
interaction takes the form 



Here L is a left-handed weak SU(2) lepton doublet; for definiteness let as assume that it 
is the first generation doublet, L = while the couplings to the second and third 

generations are suppressed by a larger scale. From the point of view of neutrino detection 
it does not matter which generation has the dominant dominant coupling to dark matter: 
the neutrino oscillations will make sure that all 3 neutrino flavors are equally populated in 
the neutrino flux from dark matter decay. 

The interaction flTT]) breaks the accidental Z2 A-parity and allows the dark matter to 
decay. It is a dimension six non-renormalizable operator, but it can be easily obtained from 
a renormalizable UV theory as follows. Let us introduce right-handed neutrinos Aj which 
have the usual Standard Model couplings. In addition, introduce a pair of fundamentals 
in the hidden sector: a scalar Q and a (Dirac) fermion x, which have masses equal to An- 
The lagrangian 



4 Dark Matter Decay 



Ant = T^AV^'^G^l.iJL + h.C.. 



(11) 



■N 



c = c 



kin 



X^QH'^x) - c^QhN^ - Q^Q\H\^ - Y^fN.HL, + h.c. 



(12) 
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contains the most general gauge-invariant renormalizable operators which can be built from 
the fields in the model. Integrating out the heavy Q, x, and Ni at one loop gives the effective 
interactions 

C^nt = -^X''a^"'G%HLL + "-^X^X'^m^ + ^G^l^G-^^H^r + h-c. (13) 

The coefficients c, C2, and C3 are of order l/167r^. The latter two interactions have negligible 
effect on phenomenology, but the first is the desired coupling of dark matter to the standard 
model via the neutrino portal. If the scale A^r is of order lO^'^ — 10^^ GeV — the typical 
scale for lepton number violating interactions — then the dark matter lifetime is of order 
10^^ seconds. Such a dark matter particle is stable on cosmological time scales and, at the 
same time, the decay rate may be large enough to allow for detection of dark matter decay 
products. 

The perturbative vertex given by (ITTl) yields two distinct partonic contributions to dark 
matter decay: first, the 1 — > 3 process X ^ gH L with a (hidden) gluon in the final state, 
and second, the 2 — > 2 process X g —>■ H L, with a (hidden) gluon in the initial state. Both 
partonic processes are part of the full gauge invariant amplitude for the decay of the dark 
matter particle, but, as we will see, the first process is the dominant contribution. The 
main contribution to dark matter decay is then the decay of the heavy adjoint A" via flTTl) to 
two-body and three-body final states with a hidden gluon plus standard model neutrinos, 
gauge and Higgs bosons. The hard gluon g from the perturbative decay and the remnants 
of the brown muck that made up the glueballino go on to hadronize into some number of 
hidden glueballs. 

We first consider the two-body decay with one neutrino and one dark gluon in the 
final state. The interaction vertex is obtained from f|TT]) by replacing the Higgs field by its 
vacuum expectation value. The partial decay width is given by 

r(A - gu) - r,, = ^ (14) 

The CP conjugate decay to anti- neutrinos has the same decay width. 

Next, consider the three-body decay to gluon, neutrino, and on-shell Higgs boson. The 
partial width for this process is 

TYl^ 1 

r(A ^ ghu) = r^-^- (3 + 44r - 36r^ - 12r^ + + 12r(2 + 3r) Inr) , (15) 

(167r)'^ A]^ 3 ^ 

where r = m\/m\ ^ 1. Again, the CP conjugate decay to anti- neutrinos occurs with the 
same rate. In the limit r — 

9 ^ 

r(A ^ ghu) = r(A ^ ghu) ^ j^jr = r3.. (16) 
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In addition, there are three-body decay modes to (i) gluon, electron, and W boson; (ii) 
gluon, neutrino and Z boson. In the hmit mw,z/'>Tix these partial widths satisfy the 
relation 

r(A ^ gZu) = r(A ^ gZu) = ^T{X ^ gW+e-) = h{X ^ gW-e+) = T,, (17) 
Ignoring gauge and Higgs masses, the dark matter lifetime is approximately 

while the ratio of two- and three-body partial widths is 



Fsb 1 m 



i^omf^Y (19) 



Varying m\ in the range 1-5 TcV, the threc-to-two-body ratio varies from 1 percent to 
30 percent. Therefore the neutrino flux from dark matter decay dominates over charged 
particle fluxes, especially when the glueballino is not much heavier than 1 TeV. 

Apart from the process where the adjoint fermion decays into a dark gluon and standard 
model fields one can also consider the annihillation of the fermion with a valence gluon. 
We argue however that this gives a negligible contribution to gluebaUino decay. First, it 
is readily apparent that the rate for the gluino to completely annihilate with its cloud of 
gluons is vanishingly small. In order for the glueballino to annihilate into the vacuum, all 
of the momentum distributed among the brown muck must be concentrated in a single 
partonic gluon which then annihilates with the heavy fermion A, or in other words, all 
of the momentum carried by the gluons must be entirely localized within the Compton 
wavelength of the heavy fermion A. The probability for this to occur is negligibly small, 
of order {A/m\)^ at most. Moreover, the probability for the operator \°'a^'^G^^ to connect 
the initial state glueballino to a single glueball at rest is suppressed. This process can be 
thought of as an annihilation of the partonic A" with a coherent sum of soft gluons inside 
the hadron. The residual soft gluons after annihilation must carry enough energy to reform 
into a stable glueball. As the glueballs have masses of a factor of a several larger than the 
strong coupling scale, an annihilation of this form is kincmatically forbidden unless the A" 
has a virtuality likewise several times larger than the strong coupling scale. Computing the 
cross-section for X"- + g ^ hu and taking pg ~ A, we can estimate 

c^rn.A(l-r)-_ (20) 
A^ 247r ^ ^ 

Strictly, as the gluon participating in the interaction is soft we need to integrate this 
expression over the field configuration of the glueballino. To obtain an estimate for the 
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annihilation rate we simply estimate the result of this integral as A. This yields an estimate 
for the rate 

mAA^(l — r)^ 



247r 



(21) 



this is suppressed relative to the two-body decay rate (fT4|) by a factor of (A/t;)^. However 
the rate for the full process should be suppressed even further, as the gluino must have 
emitted a hard gluon before participating in the annihilation in order to be sufficiently off- 
shell. This emission suppresses the annihilation process by an additional factor of a^. In 
the following, we will neglect this glueballino decay mode in comparison with the two-and 
three-body decay modes coming from the decay of the gluino alone. 



5 Interlude: Glueball Decays 

The hidden sector glueballs cannot be stable on cosmological time scales as otherwise 
they would overdose the universe. In this section we discuss the thermal relic abundance 
of glueballs and present several scenarios for their decay. We will show that there are 
multiple ways to eliminate the overabundance of glueballs, and will indicate the observable 
signals associated with each scenario. At the end of the day it will be clear that the signals 
and resulting constraints are very model-dependent, in contrast to the sharp predictions 
from two- and three-body decays via the neutrino portal discussed in Section HJ Therefore, 
when discussing signals and constraints in the following of this paper we will focus on 
the more robust signals of the standard model particles coming directly from the neutrino 
portal operator \°-G°'^^o^^ HL. The goal of this section is to point out and characterize 
a range of options for addressing the glueball problem, and demonstrate that solutions 
are possible. Readers less interested in model building and more interested in model- 
independent astrophysical signatures of the neutrino portal are encouraged to jump directly 
to the next section. 

First, we demonstrate why the glueballs must decay. Treating hidden sector confinement 
as adiabatic, the energy and entropy in the thermal gluon bath are wholly transferred to 
glueballs after confinement. If the glueballs do not decay, the present-day energy density 
in glueballs is then given in terms of the present-day photon temperature Tq by 

Pglueito) ~ ^ A , (22) 

-ju g*SM 

where g^.fl is the effective number of relativistic standard model degrees of freedom today, 
and g*,D,SM are the effective number of relativistic degrees of freedom in the dark sector 
and in the standard model when the two sectors were in thermal equilibrium. Avoiding 



13 



overclosure then requires the strong couphng scale A to be smaller than an eV, 

iN'-l)^<l. (23) 

However, the strong coupling scale A is constrained by astrophysical bounds on the dark 
matter self-interaction cross-section. Analyses of galactic cluster dynamics combine to yield 
a bound on the total dark matter self- interaction cross-section [HI [171 [IS [19] . 

— < 0.3 - 1 cmV^ ~ 1400 - 4600 GeV-^ (24) 
m 

The most stringent such astrophysical bound comes from observations of the Bullet Cluster, 
1E0657-56 [19] . In our model, the total self-interaction cross-section for the glueballino 
is geometric, (Xtot ~ vr/A^. With A ~ eV, the bound requires the mass of the glue- 
ballino to be very large, mx > 10^ GeV, which would be incompatible with thermal relic 
abundance. 

As the glueball masses are intrinsically tied to the confinement scale, there is no way to 
escape this problem. The glueballs must therefore decay. There are then three possibilities: 
the glueballs can decay within their own sector; they can decay to the standard model; or 
finally, they can decay to yet another sector. Adding light flavors to the hidden sector 
gauge theory destabilizes the glueballs to hidden sector mesons and opens an interesting 
range of possibilities for model building, which we discuss below. 

Another option, which we will pursue here, is to add additional degrees of freedom 
("mediators") which couple the hidden and visible sector, enabling glueballs to decay to 
standard model final states. Dark matter decays then yield a model-dependent amount 
of visible standard model particles from glueball decays in addition to the primary gluino 
decay, which must be taken into account. In order to avoid cosmological constraints, it is 
simplest to ask that the lifetime of at least one species of glueball be less than one second, so 
that the glueballs can safely decay without a trace before nucleosynthesis. For definiteness, 
we will consider two simple models where the mediators couple to the standard model Higgs 
boson, thereby enabling the 0"*""^ glueball to decay rapidly through the dimension 6 Higgs 
portal operator to pairs of standard model fermions. In the first model, we introduce a 
standard model singlet scalar field transforming as a fundamental under the dark gauge 
group, with the interaction 

C,r,t = y\M'\H\^ (25) 

In the second model, we introduce a pair of Dirac fermions, Kl and Pr, which are funda- 
mentals under the dark gauge group and have the standard model quantum numbers of a 
left-handed lepton doublet and a right-handed charged lepton, respectively. The interac- 
tions of the new fermions and Pr we take to be C- and P-preserving, for simplicity: 

C = /:kin + {vPrHKl + h.c.) . (26) 
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Integrating out the mediators induces the dimension six interaction 

^e// = ^i^^^^G«.G-^ (27) 

Here M is the mass scale of the mediators. The coefficient c is ?/^/12 for the scalar mediator 
model of equation (!25|) . and 2?/^/3 for the fermion mediator model of equation (p6l) . 

The 0"'""'" glueball decays through the effective interaction (1271) with branching fractions 
dictated by the couplings of the standard model Higgs. For the parameter range we con- 
sider, the O"*"^ will decay to pairs of standard model fermions through an off-shell Higgs. 
The width for this decay is given by 

= {lQ7^m\ml-ml)) ^'^'//(^o) (28) 

Here ^^^jjijn^) denotes the partial width oi h ^ ff for a standard model Higgs boson 
with mass mo, and Fq+ is the O"*""*" decay constant. We will use the lattice result for 

m , 

F(f = /(0|G;^^G'^'^'^|0++) = 6.12 m^. (29) 

If the 0"'"'*' is heavy enough to decay to bb pairs, a lifetime less than one second is readily 
accomplished with mediator mass scales in the range of tens to hundreds of TeV. Lighter 
glueballs, with masses in the range 2m,- < mQ++ < 2mb, are also viable if the mediator mass 
scales are in the range of a few to tens of TeV. Glueballs light enough to decay dominantly 
to muon pairs are interesting, as this mass range appears best compatible with obtaining 
the correct relic density for the heavy hadronized adjoint, as well as with possible visible 
signals. However, realizing a sufficiently short lifetime for the O"*""^ in this case requires 
the mediators to be so light that they lie outside the range of validity of our calculations. 
Models with light standard model-charged mediators are ruled out, but a light scalar mixing 
with the standard model Higgs is entirely possible, and seems like an interesting avenue to 
pursue. 

In the case of a scalar mediator, the only interactions induced upon integrating out the 
heavy mediators are of the form 



C^jf = G;^G'^^^ Y.J^(H^H)\ (30) 

k 

The Lorentz decomposition of the operator G'^^^G"''^'^ contains only the O"*""*". Therefore 
only the O"*"^ can decay directly via the Higgs portal to standard model particles. Other 
glueballs can decay through the O^"*" channel, e.g., 2++ (0++)*0^'^ ffO'^^- These 
three-body decay rates are highly suppressed relative to the 0+^^ decay rate due to the 
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smaller available phase space. The only glueballs which cannot decay in this manner are 
the lightest states in the P and C-odd sectors, the 0""*" and 1"*"". Without explicit P or C 
violation, no dimension 6 operator permits the 0""*" and states to decay into standard 
model fermions. 

In the case of bifundamental mediators, loops of the heavy particles also induce dimen- 
sion 8 effective interactions coupling the hidden gluons to standard model gauge bosons. 



W/ = Tri^^ Tr + F, Tr G^ (31) 



considered in detail in pT] in the context of hidden valley signals at the LHC Here F (G) 
is the field strength tensor for standard model (dark sector) gauge fields. The field strength 
tensors are contracted according to different irreducible representations of the Lorentz group 
which for the sake of brevity are not shown here. These dimension 8 operators allow all 
the hidden sector glueballs to decay, either directly to standard model gauge boson pairs 
(the O^"*", 2++, ^, 2^+) or radiatively to another glueball via the emission of a photon 
(all others). 

For a glueball which can decay either through radiative photon emission or through an 
off-shell Higgs, the situation is rather involved. The branching ratios of these two modes 
depend critically on the mass scale M of the mediators and the nonperturbative transition 
matrix elements. The ratio of the partial widths can be written parametrically as 

aadX —/ ^ (32) 



y niA ' niA 

where x = M^(l/m^ — 1/mj,) and a is the fine structure constant. The dimensionless 
phase-space function / (a, b) is given by 

I{a,b)= dxx {g{a\x-l)f^^'^ {g{h'',h''-x)f'^ . (33) 

where g{x,y; z) = {1 — x / z — y / zY — Axy / and L is the relative orbital angular momentum 
of A and B. Given the masses in Fig. [H J ^ 10~^ — 10~^. For x ~ 0.1 — 0.5, which is 
a conservative choice if the mass splittings between the mediator fields are not especially 
large, the larger suppression of dimension eight operators tends to suppress the radiative 
decays in favor of the three-body decays. However, for small enough M, the radiative 
decays can become competitive. Indeed, a simple estimate suggests that for M > 10 TeV, 
the dimension 6 operator dominates, while for mediators within the reach of the LHC 
radiative and three-body rates can be comparable. If the off-shell Higgs needs to decay 
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to muon pairs, the small muon Yukawa coupling will additionally suppress the three-body 
rates, making the radiative decays dominate. 

To summarize, glueball decays contribute b quarks or r leptons to the visible signals of 
dark matter decay, as well as possible photons (both singly and in pairs) if the mediators 
are charged under SU{2)l x U{1)y- (Decays into W and Z bosons are kinematically 
impossible.) Colored mediators would also allow the glueballs to decay into standard 
model gluon pairs, but in order to minimize the antiprotons coming from glueball decays 
we have taken the mediators to be color singlets. Moreover, as most of the mass splittings 
between glueballs arc smaller than mo, if the 0^+ is near the threshold to decay to bb pairs, 
then kinematics may prevent bb pairs from being produced in three body decays of heavier 
glueballs. A typical final visible state would then be of the form bbcc, bbr^r^ and so on. 
A similar story holds for 0+"'" decaying to tau pairs: higher-mass glueballs could yield pion 
and muon pairs in addition to the tau pairs coming from the decay of the terminal 0"*""*". 
Glueballs which decay mainly into muon or pion pairs might be obtained with additional 
model building in the Higgs sector. 

While it is possible to work out glueball branching fractions for a given mediation 
model with the help of lattice data for N — 3, it is considerably more difficult to arrive 
at the full spectrum of visible standard model particles produced in a single glueballino 
decay. The injection spectrum dNsu/dE of visible standard model particles depends on 
the number and spectra of glueballs initially produced in the glueballino decay. In order 
to make crisp predictions about indirect detection signals coming from dark matter decay, 
we therefore need to understand the process of fragmentation in pure Yang-Mills theory. 
Unfortunately, this is a situation where neither phenomenological examples from QCD nor 
data from the lattice can be of help. In the absence of light fundamentals, the color tube 
connecting hard partons cannot break. Glueball formation in a pure Yang-Mills theory 
occurs heuristically through the crossing of the tube onto itself, a process qualitatively 
different from fragmentation in QCD. 

One simple approach to estimating the relative abundances of different glueball species 
is a thermal model. In such a model, the ratio of the multiplicity of glueballs is given by 



where mj (Nj) is the mass (multiplicity) of a glueball with spin J . The effective tem- 
perature T can be taken as the center of mass energy of the colored system immediately 
following the decay of the heavy adjoint fermion, T ~ y/Am^. With mx » A, the thermal 
model predicts glueball democracy in the final state: for mx = 1 TeV and A = 1 GeV, 
the heaviest of the stable glueballs, the 0"*"", is produced two-thirds as frequently as the 




(34) 
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0+^ . Unfortunately, the applicability of the thermal model to pure Yang-Mills theory is 
unclear. In particular, excitations of the flux tube (such as can easily be produced from 
gluon showering in the final state) might enhance the production of the higher-spin states 
above the thermal estimate. Moreover, the thermal model is not well-suited to address a 
key question for indirect detection: what is the average total number of glueballs per decay? 

To arrive at an injection spectrum dNsu/dE which can be compared to indirect de- 
tection data, we need to specify not only a particular model for glueball decay but also a 
model for fragmentation in the hidden sector. While it appears possible to build interesting 
models which have the potential to avoid overproduction of antiprotons, it is clear that not 
only is there significant freedom in the detailed choice of model, but also that a reasonable 
model of pure-glue fragmentation is a necessary ingredient in quantifying the predictions 
of such models for indirect detection experiments. 

5.1 Model Building with Light Hidden Flavor 

Adding light (m^ < A) fundamental flavors to the hidden sector destabilizes the glueballs 
to hidden mesons, without destabilizing the adjoint fermion A. As A is fermionic, it must 
decay to an odd number of hidden quarks and anti-quarks. But such a final state carries 
non-zero hidden-sector baryon number. The decay of A to light flavor therefore violates 
baryon number, and is therefore forbidden in the absence of hidden sector baryon number 
violation. 

The meson masses depend on the input light quark masses as well as on the strong 
coupling scale, and can be tuned independently of the other parameters of the theory. 
Avoiding overclosure requires that the hidden mesons have sub-eV masses, giving rise to 
a large but technically natural hierarchy, ^ A ^ nix- However, the existence of such 
extremely light mesons leads to concerns that inelastic glueballino scattering \E'\1/ + 

nvr^, may exceed observational bounds on dissipative dark matter interactions. Adding such 
extremely light mesons will also dramatically increase the late-time glueballino annihilation 
cross-section, likely to a problematic extent. 

More attractive model building possibilities are realized by keeping ~ A and allowing 
the hidden sector mesons to decay to the standard model. This can be realized in multiple 
ways. One option is to introduce a Z' coupling to both hidden flavor and visible standard 
model fermions [7]. This allows at least one hidden sector meson to decay to a pair of 
standard model fermions, with a preference for heavier fermions due to helicity suppression. 
The main challenge is then obtaining a short enough lifetime for the shortest-lived hidden 
sector meson while keeping ~ A ~ m^^ at the GeV scale. For models of the form 
considered in [7], lifetimes of less than a second can be achieved for 2mj. < rriq ~ A, with 
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the mediator mass in the range of a few to tens of TeV. Decay to muon pairs can only be 
reahzed for a hght mediator and m7r„ close to the r threshold. 

Models with light flavor have one especially nice feature: they allow for more concrete 
predictions for the visible signals from dark matter decay, as with the addition of light 
flavor it is possible to adapt models of fragmentation in QCD. This then allows definite 
predictions for the spectra and multiplicities of the hidden sector hadrons produced in 
dark matter decay, which in turn allows detailed comparison to astrophysical data. While 
this is appealing, it comes at a cost: the existence of light degrees of freedom in the 
hidden sector exacerbates the uncertainties in evaluating the annihilation cross-section for 
glueballinos, and hence in the acceptable mass range for the dark matter. Nonetheless, 
there are interesting avenues for model building here. 

6 Constraints from indirect detection 

We turn to discussing the constraints on our neutrino portal scenario from various indirect 
detection experiments. The glueballino decaying via the operator of Eq. produces dom- 
inantly neutrinos, but also Higgs particles, W and Z bosons, and electrons and positrons. 
These fluxes of standard model particles produced directly from the neutrino portal opera- 
tor are relatively model- independent: the relative normalization of the two- and three-body 
decay modes depends only on the kinematics of the decay. Bounds on these decay modes 
therefore constrain the general scenario where structure in the dark sector enables dark 
matter decaying through the neutrino portal to yield a primary flux of monochromatic 
neutrinos. We compute the fluxes of neutrinos, antiprotons, positrons, and gamma rays 
observed at the Earth for standard model particles produced in the decay of the heavy dark 
matter particle, and determine the allowed parameter space of our model. 

Throughout this section we assume the dark matter distribution in our galaxy is de- 
scribed by the Einasto profile 



where r is the distance from the center of our galaxy, and we use the values a = 0.17 
and Tg = 20 kpc [22]. The normalization constant is po = 0.06GeV/cm^ which gives the 
local dark matter energy density Pq = poMi^Q) = 0.3 GeV/cm^. Here = 8.5 kpc is the 
distance of the Solar System from the galactic center. Because for decaying dark matter 
the signal depends on p£,M (rather than p|)jvf as for annihilation) there is limited sensitivity 




(35) 
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PDM{r) = po{-y 1+ -] ] (36) 



to the choice of the profile. For example, another standard profile, the NFW profile 

r/ \ \rsj , 

with parameter choices {a,/3,7} = {1,3,1}, = 20kpc, and po = 0.26GeV/cm^ 
increases the neutrino flux above that predicted using the Einasto profile by about 10 
percent. 

6.1 Neutrino bounds from Super-Kamiokande 

Dark matter decaying in our galaxy produces mostly neutrinos and antineutrinos with 
energy equal to rnxjl. This leads to a flux of cosmic neutrinos concentrated in the direction 
toward the galactic center. The interaction cross-sections of dark matter with nuclei in our 
model are negligible, leading to correspondingly negligible dark matter capture rates in the 
Earth and the Sun. Therefore only bounds on neutrino flux from the galactic center are 
constraining. To date, the best bounds on such a flux are those from the Super-Kamiokande 
collaboration [2^ [25]. 

The flux of neutrinos coming from dark matter decay is given by 

^"^'^ -(mx) = — / dsp{s), (37) 



dEydVt 4:71 mx dEy 

where dN/dE^ is the spectrum of neutrinos generated in the dark matter decay, and the 
integral of the dark matter energy density is taken along a line of sight. It is conventional 
to define the dimensionless quantity 

dJ _ f dsp{s) ^^g^ 



dn J tqPq 

normalizing the line of sight density integral to our local position and density. The flux in 
a cone of half-angle ip around the galactic center is then 

dEi, 4tt nix dE^ 

with 

r° I 

AJ(V^) = / dcosO dx p{Jrl + x"^ -2rQXCose). (40) 

icosi/) Jo 

Notice that Eq. (l39l) is the total neutrino flux, summed over all three flavors, as oscillation 
over galactic scales will equally populate all three flavors. Moreover, there is an equal in 
magnitude anti-neutrino flux. 
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Having fixed the halo profile in Eq. (!35|) . A J is a function of the half-cone angle which 
can be computed numerically. Neutrinos produced in a decay travel to the Earth without 
experiencing further interactions. The normalization of the neutrino flux at the Earth 
depends on the dark matter lifetime and mass via the ratio T2b/fnxi and this is the quantity 
that will be constrained by Super-Kamiokande. 

The published Super-Kamiokande bounds do not apply directly directly to the neutrino 
flux but rather to the upward-going muon flux. These muons are produced by charged 
current interactions of the cosmic neutrinos in the rock below the detector. In order to 
translate the neutrino and antineutrino fluxes into a muon flux we need to compute the 
conversion functions Cy^^-{E^; Ey) and Cp^^+{E^] E,^), 



dEp^ J dE^ dEy ' dE^ J dE^ dE^ 

which encapsulate the probability that a neutrino will interact in rock to produce a muon 
of energy E^ which reaches the detector. The conversion functions depend on the the muon 
range in rock R^{E^), the proton and neutron number densities in rock uni N = n,p, and 
the neutrino inelastic scattering cross section on nucleons ct^at: 



dE^ 



dfJi/p d(Ty^ 



dE^ 



dfJyp d(Tyji 



(42) 

To compute the neutrino-nucleon scattering cross sections, we use the tree-level partonic 
cross-section integrated over the CTEQ5 parton distribution functions. We perform our 
computations using the fiducial material "standard rock" , for which the nucleon densities 
are rip = n„ ^ 6.1 x 10"^^ GeV^ 

The last piece we need to evaluate Eq. ( l42l) is the muon range i?^(£'^), the distance a 
muon with the initial energy E^ can travel before its energy is radiated away so as to fall 
below the Super-Kamiokande threshold Eq = 1.6 GeV. The muon range can be modeled 
analytically using the "continuous slowing down approximation", with parameters found 
in ^26|. We approximate the muon range R^{E^) as 



R^{E^) ^ lO^cm In 



E^ + e 



(43) 



Here e is the critical energy, where ionization and radiation loss are equal; we take e 
to be 693 GeV in standard rock [26] (we assume that cross sections in standard rock are 
comparable to cross sections in the rock below Super-Kamiokande). The continuous slowing 
down approximation becomes poor at high energies, where stochastic radiative processes 
dominate the energy loss. More accurate Monte Carlo modeling of muon propagation 
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Figure 4: The maximal V2h/'mx as a function of mx allowed by the Super-Kamiokande 
bounds. 

at high energies reduces the flux of upward-going neutrino-induced muons at neutrino 
detectors for a given incident neutrino flux relative to the flux obtained using the continuous 
slowing down approximation [27J. Comparing the effective muon range given by (H3!) to 
the effective muon range computed in [27] for muons with energy > TeV, we find our 
muon ranges are overestimated by a factor of less than order unity over most of the energy 
range. The bounds we obtain from Super-Kamiokande's reported upward through-going 
muon counts are therefore conservative. 

Finally, to obtain the total muon flux in the Super-Kamiokande detector we have to 
integrate over all muon energies between Eq and vrL\j1. At the end of the day the muon 
flux in Super-Kamiokande is given by 

The factor 1/3 arises because 2/3 of the produced neutrinos oscillate into electron and tau 
neutrinos to which Super-Kamiokande is far less sensitive. 

The muon flux induced by neutrinos from the galactic center predicted by our model is 
dominated by the monoenergetic neutrinos from the dominant two-body decay. At higher 
dark matter masses, rrL\ ~ 5 TeV, the additional neutrinos from the sub-dominant three- 
body modes begin to make a noticeable contribution to the total muon flux. We computed 
the total muon flux predicted from our model and confronted it with the Super-Kamiokande 
bounds from |2l] . The most stringent bounds are obtained for cones of half-angle 30 degrees. 
In Fig. mthe resulting bound on V^hl'^x is plotted as a function of m\. The dark matter 
mass is varied in the range 1 — 5 TeV, as motivated by the estimate of thermal abundance 
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in the previous section. For m\ = 1 TeV the bound is T2b/ni\ < 3.2 x 10^^^, corresponding 
to the hfetime r 9.9 x 10^^ seconds. For m\ = 5 TeV the bound is T2b/^\ ~ 1-6 x 10~^^, 
corresponding to the hfetime ~ 1.9 x 10^^ seconds. 

The sub-dominant three-body decays lead to additional visible signals, in particular, 
to antiprotons, hard leptons (which we take here as e^), and photons. In the following 
subsections we investigate the predicted antiproton, positron and gamma ray fluxes arising 
from the three-body decay modes to determine whether these subleading decays impose 
more stringent bounds on the lifetime of our dark matter particle than the two-body decays. 

6.2 Antiproton bounds from PAMELA 

An irreducible source of antiprotons in our model are the decays of the W , Z, and Higgs 
bosons which are produced by the three-body decays of our dark matter particle. On 
average, these channels lead to one energetic antiproton per decay. The branching ratio 
for three-body decays varies from 4 percent of the total width for m\ = 1 TeV to as much 
as 55 percent for mx = 5 TeV. Thus we expect that the antiproton bounds on the dark 
matter lifetime can become quite stringent in some regions of the parameter space. 

The production rate of antiprotons per unit energy due to the three-body decays is 
given by 

Q,iE, r) = 2-poM{r) [-^ + ^ + J • (45) 

We determined the antiproton spectra dN.^p/dE by first generating W, Z, and Higgs 
decays at rest using Pythia 6.420, and then boosting the momenta of antiprotons to the 
galactic rest frame. The energy distribution of the parent gauge and Higgs bosons is 
given by the differential decay width, which is the same for all three particles in the limit 
mx > mh,w,z- 

1 dT^b ^'^.E(mx-2E){mx-E). (46) 



Tsb dE ml 

Throughout, the Higgs is assumed to have standard model properties and the mass equal 
115 GeV. 

In order to translate the production rate into the flux of antiprotons at the Earth it 
is necessary to model the propagation of antiprotons in our galaxy [28]. The antiproton 
number density per unit energy Np satisfies a diffusion equation 

K{E)V^Np{E,r) - d,{V,sgn{z)Np{E,r)) - 2h6{z)Tann{E)Np{E,f) + Qp{E,r) = (47) 

in which the production rate (1451) appears as the source term. Energy loss for protons is 
negligible and the corresponding term in the diffusion equation is dropped. The diffusion 
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Figure 5: Antiproton flux from tlie tliree-body decay of the gauginos (solid red) compared 
to the flux measured by PAMELA (dashed) for dark matter masses mx = 1, 2, and 5 
TeV. In each case the glueballino decay rate is set to the maximal value allowed by the 
Super-Kamiokande bounds. 

coefficient K[E) is parameterized as K[E) = Kq{E/ GeV)*^. The parameters Kq and S 
are not known from first principles; instead, they are fitted, along with the magnitude of 
the convective wind Vc and the size of the diffusion zone, so as to correctly describe the 
observed flux of cosmic rays. The third term describes the disappearance of antiprotons due 
to annihilations with protons in the galactic disk. The half-height of the galactic disk is set 
hj h = 0.1 kpc, and the annihilation width is given by the number densities of hydrogen and 
helium multiplied by the proton-antiproton annihilation cross-section, TanniE) = [uh + 
A^^^nHe)(^pp{E). Expressions for nH-,nHe and the proton-antiproton annihilation cross- 
section (Tpp{E) can be found in [2S]. The diffusion equation Eq. fH7|) is solved subject to 
the boundary condition that the antiproton density Np vanish at the surface of the diffusion 
region, taken to be a cylinder of radius Rmax = 20 kpc and half- height Lmax- 

In the following we focus on the MED propagation model [30], in which case the diffusion 
equation parameters take the following values: 

K{E) = 3.5 X 10"^°kpcVs; 5 = 0.7; V^=1.2x lO^cm/s; L^nax = 4 kpc. (48) 

Once we know Np, the flux of antiprotons observed at the Earth is given by = 
^iVp(£', -Tq); in the interesting energy range the antiprotons are highly relativistic and 
V ^ 1. 

The propagation equation can be solved semi-analytically |3T]. Using that solution, the 
antiproton flux in our model is given by 

1 Tsb f dNz-.p dNh^p dNw-*p\ , . 

The propagation effects are encoded in the function R{E) which can be computed numer- 
ically [an [30]. 
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With this solution in hand we are ready to estimate the antiproton constraints on the 
hfetime of our dark matter particle. We pick three representative values for the dark matter 
mass: rrix = 1,2,5 TeV. For the sake of comparison with the neutrino bounds, for each 
mass we choose the maximum T2b/fn\ allowed by Super-Kamiokande as plotted in Fig. HI 
Given r2f, and m^, the partial width for three-body decays Fs;, is completely determined 
by equation fll9p . The mass fixes the spectra dN/dE, while the partial width fixes the 
normalization of the antiproton flux in Eq. ( l49l) . Our results are plotted in Fig. [5] and 
compared to the antiproton flux measured by PAMELA. It is clear that antiprotons, despite 
being produced by subdominant decay channels, can put important constraints on the dark 
matter lifetime. For mx = 5 TeV, the predicted antiproton flux is 5 times larger than 
observed, which would imply the stringent bound T2b/m\ ^ 3 x 10"^^ corresponding to the 
lifetime ~ 10^^ sec. For m\ = 1 TeV, the predicted antiproton flux is only marginally larger 
than observed, and is slightly more constraining than the Super-Kamiokande bound. One 
should however keep in mind that the propagation of antiprotons in our galaxy suffers from 
huge uncertainties, and the predicted antiproton flux changes dramatically with different 
assumptions about the propagation model. For example, choosing the MIN parameters for 
the propagation model [30] . 

Ko = 5.1 X lO^^^kpcVs; 5 = 0.85; K = 1-35 x lO^cm/s; L^nax = 1 kpc, (50) 

reduces the antiproton flux by a factor of five. This would make the antiproton bounds 
comparable to the neutrino ones for nix = 5 TeV, and irrelevant for mx = 1 TeV. (It is of 
course also possible to choose propagation parameters which magnify the antiproton flux.) 
Thus, the antiproton bounds should be taken with a grain of salt. In the following, we 
continue using the more robust Super-Kamiokande constraints as the reference point for 
estimating the indirect signals of dark matter. 

6.3 Other constraints: photons, positrons 

Our decaying dark matter particle contributes also to the cosmic flux of positrons and 
photons. Hard positrons are produced together with W bosons in three-body decays and, 
to a lesser extent, by subsequent decays of W, Z, and Higgs. Gamma ray photons are 
produced via decays of W, Z, and Higgs, from final state radiation off charged particles, 
and through inverse Compton scattering of the hard positrons and electrons on starlight. 
Below we compute the fluxes of hard positrons and photons from the three-body decay 
modes and argue that the corresponding bounds are far weaker than those from neutrinos 
and antiprotons. Therefore, the PAMELA positron measurements and the gamma ray 
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measurements from FERMI and HESS do not impose new constraints on the lifetime of 
our dark matter particle. 

The positron flux at the Earth due to the positrons from three-body decays is given by 

m 



Here, the energy loss coefficient b{E) gives the rate of energy loss due to scattering off 
of cosmic photons and magnetic fields, and is approximately given by b{E) = teE"^/ GeV 
where te = lO^^s. The function I{E; E'), sometimes called the halo function, results from 
semi-analytically solving the diffusion equation for positrons [32]. It is analogous to the 
function R{E) for antiprotons, and similarly depends on the choice of propagation model. 
Finally, dNe+ / dE is the positron spectrum from three-body decays, given by 



^dV^ ^ 32 / ^ _ 2^ 
V-xh dE mi \ 3 



E^mx--E]. (52) 

To obtain the positron fraction one has to divide the flux in Eq. (|5T|) by the total flux 
of electrons and positrons, for which we take 

(i$e± 0.012 , 2 r. .roN 

= l/(cm2 sec GeV) (53) 

as measured by FERMI [3]. As an example, in Fig. [6] we plot the contribution to the 
positron fraction from the three-body decays for rrix = 2 TeV, using the MED propagation 
parameters. We took the decay rate T^h/mx to be the maximum one allowed by the 
neutrino bounds from Super-Kamiokande. Recall that for these parameters antiprotons 
are overproduced by a factor of 2-3. On the other hand, the positron fraction for the same 
propagation parameters is a factor of 2-4 below that measured by PAMELA. We conclude 
that positrons do not impose any new constraints on the lifetime of our dark matter particle. 
At the same time, we cannot explain the anomalous positron fraction via the positrons from 
three-body decays; another primary source of positrons is needed to explain the PAMELA 
positron measurement. The new source could be of purely astrophysical origin [53], but it 
could be glueball (or meson) decays in our model. 

Similarly, the contribution to the gamma ray spectrum from the three-body decays is 
well below the current observational bounds. As an example, we calculate the contribution 
to the diffuse gamma ray spectrum at intermediate galactic latitudes. Recently, FERMI 
published the gamma ray spectrum between 0.1 and 10 GeV from between 10 and 20 
degrees galactic latitude [3l|. The contribution of three-body decays to that spectrum is 
given by 

d<!>^ _ TsbdN^^ 1 



,^ 2^— ^J,(10°,20°) I,{^,,^,) = — d<p dl rfcos^pDAf(r) (54) 
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Figure 6: Contributions from three-body decays to positron and gamma ray fluxes for 
rrix = 2 TeV. Left: the positron fraction (sohd red) compared to the PAMELA measurement 
(dashed). Right: the diffuse gamma ray flux 10° < 6 < 20°, < / < 360° (solid red) 
compared to the FERMI measurement (dashed). 

where r = ^/tq + P — 2rl cos 9 cos 0. Using Pythia, we determined the photon spectrum 
dN^/dE from W , and Higgs decays, where the parent bosons originate from the three- 
body decays of dark matter. The results for m\ = 2 TeV are plotted in Fig. El It is clear 
that the contribution of the three-body decays is more than an order of magnitude smaller 
than the flux measured by FERMI. Similarly, the contribution to the gamma ray flux 
from the galactic center is far below the flux measured by FERMI and HESS [35]. Thus, 
the strongest constraints on the dark matter lifetime in our model are the galactic center 
neutrino measurement by Super-Kamiokande and the cosmic antiproton flux measurement 
by PAMELA. 



7 Discovery at Neutrino Telescopes 

Our model predicts a monochromatic flux of high-energy {Eu ~ TeV) neutrinos and an- 
tineutrinos from the galactic center, 

where the line-of-sight integral over the dark matter density proflle A J/^ (■?/') is given in Eq. 
( l39l) . The normalization of the flux is set by the dark matter lifetime, r2b- The flux is 
equally populated by all three species of neutrinos; the flux of muon neutrinos to which 
experiments are most sensitive is one third of Eq. (l55l) . The Super-Kamiokande results place 
constraints on the dark matter lifetime T2b/m\i as shown in flgure Fig. |H The dark matter 
lifetime is also constrained by the antiprotons produced in subdominant three-body decay 
modes. While the limits on cosmic antiprotons can be comparable to or more stringent 
than the bounds from Super-Kamiokande, depending on the choice of model for charged 
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particle propagation in the galaxy, to discuss signals at neutrino telescopes we will continue 
to use the limits on the dark matter lifetime set by Super-Kamiokande. In this section we 
discuss the sensitivity of the ANTARES and IceCube neutrino telescopes to the primary 
neutrino signal. Both experiments have the potential to discover this signal. 

7.1 ANTARES 

ANTARES is a high-energy neutrino telescope located in the Mediterranean sea. The 
flux of muon neutrinos from the galactic center produces upward-going muons via charged 
current interactions in the rock and water below the detector, which ANTARES observes 
via the Cherenkov light produced in sea water by the muons. The effective detector area 
Af.ff{E) for muon neutrinos can be found e.g. in [36| and at 1 TeV is given approximately 
by 50cm^. The energy resolution near 1 TeV is expected to be a factor of 2-3 [37], while 
the angular resolution is better than a degree and limited mostly by the kinematics of the 
charged current interactions [36] . 

The number of signal counts per second from the cone of half-angle ip around the galactic 
center is given by 

and the neutrino energy is mx/2. 

The irreducible background for the dark matter decay signal comes from atmospheric 
muon neutrinos. The combined flux of atmospheric muon neutrinos and anti-neutrinos at 
1 TeV is, averaged over angle [38] , 



atm 



dEdn 



= 2.7 X 10"''( GeVcm's si)-\ (57) 

ITeV 



The atmospheric neutrino spectrum falls as in the energy range of interest, 100 GeV < 
£' < 10 TeV. In the energy range 500 GeV < E < 5 TeV where we expect hard primary 
neutrinos from decaying dark matter, neutrino telescopes have poor energy resolution. To 
differentiate signal from background we therefore rely on first, absolute rate, and second, 
angular variation in the signal. The neutrinos coming from dark matter decay are concen- 
trated at the galactic center, while the spatial variation of the atmospheric neutrino fluxes 
is well-known and is at most a factor of 2.5 between zenith and horizon [38]. The number 
of counts per second from atmospheric neutrinos in a given energy bin Emin, E^ax is given 
by 

^^.r,EA,AE)'^ (58, 



28 



ANTARES Mnm = ITeV ANTARES Mqm = CTeV ANTARES Mnm = 5TeV 
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Figure 7: Neutrino counts per year in Antares from the signal (solid red) and the back- 
ground of atmospheric neutrinos (dashed blue) as a function of the halfangle ip around the 
galactic center, for mx =1, 2, and 5 TeV. The normalization of the flux corresponds to the 
maximal flux allowed by Super-Kamiokande for a given dark matter mass. 

We approximate the energy resolution of ANTARES by choosing an energy bin rrix/Q < 
E < 3mx/2, corresponding to a factor of three in energy resolution centered on the energy 
of the primary neutrino signal. In Fig. [7] we plot the number of signal and background 
counts per year in this bin for cones around the galactic center as a function of the cone 
half-angle. We show plots for mx = 1, 2, and 5 TeV. To generate these plots, we assume 
that the detector effective area ^e// already incorporates factors accounting for detection 
efficiency and the different cross sections of neutrinos and anti-neutrinos. We use the 
atmospheric z/^ + Up, flux d^atm/dEdQ computed in [38] and averaged over the zenith 
angle. For the signal, we show expected counts for the maximal decay rate allowed by the 
Super-Kamiokande bounds. For that signal rate, the flve-degree cone corresponds to on the 
order of 5 signal counts and a signal-to-background ratio in the range of 1-3 (for 1-5 TeV 
dark matter). Thus, in the best case scenario, a clear signal can be detected already in 
one year of ANTARES data. If no signal is detected, ANTARES can signiflcantly improve 
the bounds on the dark matter lifetime in our scenario. The theoretical uncertainty in the 
atmospheric neutrino flux is of order 25 percent at 1 TeV [38] . From that we estimate that 
one year of ANTARES data can improve the bounds on by a factor of 4 to 8, depending 
on the dark matter mass. 

7.2 IceCube 

While IceCube, located at the South Pole, has better exposure to the Northern Hemisphere, 
the extension DeepCore, planned for 2010, will allow IceCube to observe the Southern 
Hemisphere sky [32] • IceCube's new sensitivity to the galactic center opens many excit- 
ing possibihties for dark matter studies |5j. DeepCore adds additional instrumentation 
in the clearest and best-shielded portion of the IceCube detector, and enables study of 
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downward-going neutrinos by using the remainder of the IceCube detector to veto cosmic 
muons. Cosmic downward-going neutrino signals then consist of neutrino-initiated events 
with initial interaction vertices located within the fiducial volume of DeepCore. 

Muon neutrinos are primarily detected in DeepCore as they are in Super-Kamiokande 
and ANTARES, as a track from the muon produced in a charged- current interaction of the 
neutrino with matter. In addition to these track-like events, the neutral-current interactions 
of neutrinos with all flavor, as well as the charged- current interactions of electron and 
tau neutrinos, deposit energy in a localized region and are reconstructed as cascade-like 
events. Cosmic neutrinos, and in particular neutrinos from dark matter decay, populate all 
three flavors equally. Usually at neutrino telescopes, the effective volume probed by muon 
neutrinos is orders of magnitude larger than the effective volume seen by electron and tau 
neutrinos, due to the large muon penetration depth. However, when using DeepCore to 
observe downward-going neutrinos, the volume available for muon neutrino events is equal 
to the volume available for electron and tau neutrino events, putting observable signal 
events from all three neutrino flavors on an equal footing. While angular resolution for 
cascade events is less than angular resolution for tracks [32] , the high rate of cascade events 
predicted by our model suggests another potential handle for discovery. 

The fiducial volume of DeepCore is estimated to be around 20 megatons [521 SD] , that 
is Vp ~ 0.02 km^. Detection efficiencies for neutrino-initiated events are estimated to be 
20%-30% [33]. While the angular resolution that can be achieved for DeepCore-initiating 
events is still to be determined, a reasonable goal is 10 degrees for tracks, possibly improving 
to 5 degrees, and 40 degrees for cascades [H]. We will estimate the energy resolution to 
be Tf = 0.5 in logE for tracks and Vf. = 0.3 in logE for cascades [Hj. The number of 
neutrino-initiated muons in DeepCore from a given incident neutrino flux is 

^= [dEj dE.a,V,Tn, ( ^ + . (59) 

dtdn J ^ \dE^dn dEf, dE^dVL dE^ J ^ ' 

The sum runs over nucleons, N = n, p. We take ice to have a density 0.92 g/cm^, corre- 
sponding to the nucleon number densities rin ~ 1.9 x 10^^^ GeV^, Up ^ 2.4 x 10^^^ GeV^. 
We include a factor ap = 0.25 to account for the detection efficiency. The integral over 
the neutrino energy is trivial for the signal; for the background we set the lower integration 
limits at E,^ = e'^*'''mx/2 for tracks and cascades respectively. 

In Fig. [8] we plot the expected number of DeepCore-initiating track-like events originat- 
ing within a cone of given half-angle around the galactic center. For the maximal signal flux 
allowed by the Super-Kamiokande constraints, one expects on the order of 25 (5) counts 
per year inside a cone of half-angle 30° for a dark matter mass of 1 (5) TeV. The dominant 
background from muon atmospheric neutrinos, also plotted in Fig. [HI gives roughly the 
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Figure 8: Track-like events per year in DeepCore due to the signal (solid red) and the 
background of atmospheric muon neutrinos (dashed blue) as a function of the half-angle 
if) around the galactic center, for m\ = 1, 2, and 5 TeV. The normalization of the flux 
corresponds to the maximal flux allowed by Super-Kamiokande for a given dark matter 
mass. 
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Figure 9: Cascade-like events per year in DeepCore due to the signal (solid red) and the 
background of atmospheric neutrinos (dashed blue) as a function of the halfangle ip around 
the galactic center, for nix =1? 2, and 5 TeV. The normalization of the flux corresponds to 
the maximal flux allowed by Super-Kamiokande for a given dark matter mass. 

same number of track-like events as the signal. However, it may be possible to reduce as 
much as 99 percent of the atmospheric neutrino background by vetoing on the associated 
collinear muon [12]. If that proves to be the case, the background could be almost com- 
pletely suppressed and the discovery prospects for the signal are greatly enhanced. The 
sensitivity of IceCube would only be limited by statistics, and could reach a lifetime of 10^^ 
seconds after several years of operation, two orders of magnitude better than the current 
bounds. 

We also compute the number of cascade-like events initiated by the signal and compare 
to events from the atmospheric neutrino background As the atmospheric electron neutrino 
flux is approximately 5 percent of the atmospheric muon neutrino flux in the energy range 
of interest, the background is dominated by the neutral-current interactions of atmospheric 
muon neutrinos. Charged-current interactions of atmospheric electron neutrinos contribute 



31 



approximately an eighth of the total background event rate. For cascade events, rate alone 
can serve to distinguish signal from background, as the signal to background ratio can 
be much larger than one for decay rates near the Super-Kamiokande bounds. In Fig. [9] 
we plot signal and background rates for cascade events for mx = 1, 2, and 5 TeV, for the 
maximal signal flux allowed by the Super-Kamiokande bounds. Again, it may be possible to 
nearly eliminate the background due to atmospheric muon neutrinos, which would reduce 
the background by nearly 90 percent. Thus cascade-like events, even more than track- 
like events, offer powerful handles on neutrino signals of this nature. More importantly, 
observing an excess of both tracks and cascades would be an important confirmation of the 
galactic origin of the neutrino flux. 

8 Conclusions 

We presented a scenario where decay of hidden sector dark matter can give a large flux 
of neutrinos from the galactic center, dominating over signals in more standard detection 
channels. Fluxes of the more readily observed e^, antiprotons, and photons are sublead- 
ing. Observational bounds on the flux of cosmic antiprotons are stringent, and despite 
the relatively small production rate of antiprotons in this class of model, these observa- 
tions can become constraining if the lifetime of the dark matter particle is of order 10^^ 
seconds. However, the large astrophysical uncertainties in the galactic propagation model 
for antiprotons make it difficult to draw any robust conclusions from the antiproton spec- 
trum. Moreover, such a dark matter particle naturally interacts very weakly with nucleons, 
giving null predictions for direct detection experiments. In these circumstances, neutrino 
telescopes become the most sensitive probe of dark matter. 

In the example model we present here, the dark matter candidate is a fermion be- 
longing to a hidden sector which is neutral under all standard model symmetries. Such 
a particle naturally couples to the standard model via the neutrino portal, which is the 
lowest-dimensional fermionic gauge singlet operator. The coupling through the neutrino 
portal enables dark matter to decay to the standard model. The existence of lighter states 
in the dark sector allows the decay channel with a single standard model neutrino to dom- 
inate the final state. This decay mode results in a monoenergetic fiux of cosmic neutrinos 
that can be probed by neutrino telescopes. While we have studied a specific realization of 
that scenario, where the hidden sector consists of a confining SU (N) gauge group together 
with a heavy fermion in the adjoint representation, obviously similar phenomenology can 
be obtained in a wider class of models. 

Currently, the coupling via the neutrino portal is most robustly constrained by Super- 
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Kamiokande, which restricts the dark matter hfetime to be larger than 10^^ seconds. The 
PAMELA measurement of the antiproton flux yields roughly comparable constraints due to 
the antiprotons arising from the subleading dark matter decay channels containing gauge 
and Higgs bosons. These bounds leave a lot of room for discovery at the current gen- 
eration of neutrino telescopes. The good angular resolution of both ANTARES and Ice- 
Cube/DeepCore ensure that the number of signal counts from the direction of the galactic 
center may be larger than or comparable to the expected background from atmospheric 
neutrinos. In the case of ANTARES, the sensitivity to the dark matter decay rate exceeds 
that of Super-Kamiokande by a factor of several. Furthermore, IceCube/DeepCore has 
the additional prospect of rejecting the vast majority of atmospheric neutrino events by 
observing the associated muon (which is impossible at ANTARES). This implies increased 
sensitivity by another order of magnitude with respect to ANTARES. 

Our analysis of the sensitivity of ANTARES and IceCube/DeepCore should be consid- 
ered only as suggestive of the experiments' ultimate capabilities; the operational parameters 
of these experiments are subject to ongoing research by the collaborations. It is conceivable 
that designing search algorithms targeted for a monochromatic neutrino signal can further 
increase the sensitivity to our scenario. 
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